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The  direct  methanol  fuel  cell  (DMFC)  has  a  significant  potential  in  consumer  electronics  and  in  backup 
and  portable  power.  Its  progress  is,  however,  hindered,  in  part,  by  a  lack  of  an  adequate  fundamental 
understanding  of  the  effect  of  various  operating  and  design  variables  on  its  performance.  While  detailed 
computational  models  are  available,  an  analytical  model  is  attractive  for  ease  of  comprehension  and 
ready  utility.  Therefore,  we  have  developed  a  comprehensive  yet  tractable  one-dimensional,  isothermal, 
explicit  analytical  model  based  largely  on  a  priori  parameters,  and  in  terms  of  quantities  with  tangible 
meaning. 

The  model  correctly  predicts  the  extent  of  methanol  crossover  and  its  effect  on  open-circuit  voltage 
(OCV)  as  well  as  on  polarization  of  the  anode,  cathode,  and  the  fuel  cell.  It  also  accurately  describes 
the  effects  of  methanol  feed  concentration  and  temperature  on  DMFC  performance.  It  aptly  predicts 
the  significant  power  losses  from  the  large  anode  and  cathode  overpotentials  as  well  as  from  methanol 
crossover,  and  the  resulting  low  DMFC  efficiency,  except  over  a  narrow  range  of  operating  conditions. 
The  insightful  model  can  be  used,  e.g.,  in  real-time  control  of  DMFC  to  operate  in  the  narrow  region  of 
high  efficiency  and  power  density. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  direct  methanol  fuel  cell  (DMFC)  [1]  is  finding  increasing 
application  in  portable,  consumer  electronics,  and  backup  power. 
A  key  appeal  of  the  DMFC  is  its  ability  to  use  a  high  energy  den¬ 
sity,  commonly  available,  and  an  inexpensive  liquid  fuel,  namely 
methanol.  In  comparison,  hydrogen  used  by  the  conventional 
polymer-electrolyte  membrane  (PEM)  fuel  cell  has  distribution  and 
storage  challenges.  The  main  shortcomings  of  the  DMFC  are  low 
power  density  coupled  with  an  order  of  magnitude  larger  amount 
of  precious  metal  catalyst  usage. 

The  conventional  DFMC  shares  the  polymer-electrolyte  mem¬ 
brane  as  well  as  an  intricate  nanostructure  with  the  PEM  fuel  cell 
[2].  Although  great  advances  have  been  made  in  recent  years  in 
its  design,  operation,  and  performance,  further  progress  depends 


Abbreviations:  ACL,  anode  catalyst  layer;  ADL,  anode  diffusion  layer;  CL,  cata¬ 
lyst  layer;  CCL,  cathode  catalyst  layer;  CDL,  cathode  diffusion  layer;  DMFC,  direct 
methanol  fuel  cell;  ECSA,  active  electrocatalyst  surface  area;  GDL,  gas-diffusion 
layer;  MEA,  membrane  electrode  assembly;  MOR,  methanol  oxidation  reaction; 
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in  part  on  a  better  fundamental  understanding  of  its  functioning. 
Consequently,  many  theoretical  models  have  been  developed  over 
the  past  decade,  the  majority  belonging  to  two  categories:  (1 )  com¬ 
putational  models  [3],  which  involve  few  assumptions  but  provide 
only  numerical  results,  and  (2)  semi-empirical  models  [4].  Both 
these  approaches,  naturally,  provide  limited  insights.  On  the  other 
hand,  a  tractable  analytical  model,  developed  from  fundamental 
considerations  coupled  with  appropriate  simplifying  assumptions, 
has  the  advantages  of  transparency  and  simplicity  of  usage.  This  is 
our  objective  here. 

The  key  technological  issues  in  DMFC  are  well  known  [1].  For 
one,  methanol  oxidation  suffers  from  sluggish  kinetics  because 
the  CO  produced  during  the  methanol  oxidation  reaction  (MOR) 
blocks  catalyst  sites  [5].  Another  issue  is  that  the  membrane  does 
not  effectively  keep  the  methanol  and  oxygen  apart,  because  of  its 
reliance  on  water  for  proton  transport.  Of  course,  where  there  is 
water,  there  is  methanol.  The  resulting  methanol  crossover  creates 
a  mixed  potential  at  the  cathode  and  a  low  efficiency  [6]. 

As  a  result,  methanol  feed  concentration  has  a  profound  effect  on 
the  performance  and  efficiency  of  a  DMFC;  higher  concentrations 
promoting  CO  poisoning  of  anode  electrocatalyst  and  methanol 
crossover,  but  low  concentrations  providing  low  power  densities 
[7-9].  Further,  operating  temperature  can  significantly  influence 
electrode  activation  and  Ohmic  resistance  [9]  as  well  as  methanol 
crossover  [10].  PEM  thickness  increases  Ohmic  resistance,  but 
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Nomenclature 

A  MEA  geometric  area,  cm2 

ci<a  concentration  of  species  z  in  layer  a,  mol  cm-3 

Cj  o  concentration  of  species  z  at  equilibrium,  mol  cm-3 

ci  b  concentration  of  species  z  in  bulk  phase,  mol  cm-3 

c,  in  concentration  of  species  z  in  feed,  mol  cm-3 

Ci  ref  reference  concentration  of  species  i,  mol  cm-3 
d  channel  depth,  cm 

dMA  anode  metal  catalyst  nanoparticle  diameter,  cm 

dM,c  anode  metal  catalyst  nanoparticle  diameter,  cm 

Djj  binary  molecular  diffusion  coefficient  between 

species  z  and  j,  cm2  s-1 

Deia  effective  diffusion  coefficient  of  species  z  in  layer  a , 
cm2  s-1 

DMeA  diffusion  of  methanol  in  the  anode  gas  diffusion 
layer 

f^Me,w  diffusion  of  methanol  in  water 
Doe  diffusion  of  oxygen  in  the  cathode  gas  diffusion  layer 
D0(W  diffusion  of  oxygen  in  water 
Ea anode  effective  activation  energy  for  methanol  oxi¬ 
dation  reaction,  J  mol-1 

£C  (p0  cathode  effective  activation  energy  for  oxygen 
reduction  reaction,  J  mol-1 

E ^  activation  energy  for  viscosity  of  water,  J  mol-1 

F  Faraday’s  constant,  96,487  C  equiv.-1 

z  fuel  cell  current  density  (A  cm-2  of  geometric  elec¬ 

trode  area) 

iA  o  anode  exchange  current  density,  A  cm-2 
z*  0  rejr  anode  reference  exchange  current  density,  anode, 
reference,  A  cm-2  active  metal  catalyst  area 
iA  i  anode  limiting  exchange  current  density,  A  cm-2 
z’c, o  cathode  exchange  current  density,  A  cm-2 

z*  0  rejr  cathode  reference  exchange  current  density,  A  cm-2 

active  metal  catalyst  area 

z'c,l  cathode  limiting  exchange  current  density,  A  cm-2 
ix, o2  crossover  current  density  of  oxygen,  A  cm-2 

ix, Me  crossover  current  density  of  methanol,  A  cm-2 

ix,Me,L  limiting  crossover  current  density  of  methanol, 
A  cm-2 

ip  electrode  p  current  density,  A  cm-2  geometric  elec¬ 

trode  area 

i*p  current  density,  A  cm-2  metal  catalyst  surface 

ip  Q  exchange  current  density,  A  cm-2  geometric  elec¬ 
trode  area 

z*  o  exchange  current  density,  A  cm-2  metal  catalyst 
surface 

KMe  adsorption  equilibrium  constant 

I<a  equilibrium  constant  for  acid  dissociation 
ki  bulk-phase  mass  transfer  coefficient  for  species  z, 

cms-1 

k*p  rate  constant  for  electrode  reaction 
/<*  rate  constant  corresponding  to  equilibrium  elec¬ 

trode  potential 
L  channel  length,  cm 

La  thickness  of  layer  a ,  cm 

mMtA  anode  catalyst  loading,  g  cm-2 

mM,c  cathode  catalyst  loading,  g  cm-2 

n  number  of  channels  in  parallel 

Ni<z  flux  of  species  z  in  the  membrane  along  the  z  direc¬ 
tion 

Pit(X  permeance  of  species  z  in  layer  a,  mol  cm-1 
Pea  effective  permeance  of  species  z  in  layer  a,  including 

bulk-resistance,  mol  cm-1 


Pi  partial  pressure  of  species  z,  atm 

pi  ref  reference  partial  pressure  of  species  z,  atm 

q  heat  dissipation  flux,  W  cm-2 

qc o2  volume  fraction  of  C02  bubbles  inside  ADL  pore  vol¬ 

ume 

qw  volume  fraction  of  water  inside  CDL  pore  volume 

P  power  density,  W  cm-2 

R  gas  constant 

Re  Reynolds  number 

R/  interfacial  resistance 

rp  reaction  rate,  mol  cm-2  geometric  electrode  area 

r*  reaction  rate,  mol  cm-2  metal  catalyst  surface 

rp; o  exchange  reaction  rate,  mol  cm-2  geometric  elec¬ 

trode  area 

r*0  exchange  reaction  rate,  mol  cm-2  metal  catalyst 
surface 

Sc  Schmidt  number 

Sh  Sherwood  number 

T  temperature,  K 

T ref  reference  temperature,  298  K 

V  cell  voltage,  V 

V0  thermodynamic  voltage,  V 

Vmax  maximum  theoretical  cell  voltage  determined  from 
AH,  V 

VA  volumetric  flow  rate  of  anode  feed,  cm3  s-1 
Vc  volumetric  flow  rate  of  cathode  feed,  cm3  s-1 
v  fluid  velocity  in  flow  channels,  cm  s-1 

X/  mole  fraction  of  species  z 

Zj  charge  number  of  species  z 

Greek  symbols 
a  layer  a 

aA  transfer  coefficient,  anode 

ac  transfer  coefficient,  cathode 

dtp  transfer  coefficient  of  the  RDS  in  electrode  reaction 
P 

(5  degree  of  acid  dissociation  in  Nation 

yMA  roughness  factor,  ratio  of  active  anode  catalyst  sur¬ 
face  area  to  MEA  area 

yM>c  roughness  factor,  ratio  of  active  cathode  catalyst 

surface  area  to  MEA  area 

S  ratio  of  mutual  to  matrix  effective  diffusion  coeffi¬ 

cients 

AG  Gibbs  free  energy  change  J  mol-1 

AH  enthalpy  change  J  mol-1 

s  fuel  cell  efficiency 

ea  void  fraction  of  layer  a 

rjA  anode  overpotential,  V 

Pb  PEM  overpotential,  V 

Pc  cathode  overpotential,  V 

?7/  interfacial  overpotential,  V 

pp  overpotential  of  electrode  reaction  p  =  O  -  Op  0»  V 
px<A  Anodic  overpotential  at  cathode,  V 
px,c  Cathodic  overpotential  at  anode,  V 

0C o-s  fraction  of  catalyst  sites  occupied  by  CO 

Ocosjef  reference  fraction  of  catalyst  sites  occupied  by  CO 
Ki<a  partition  coefficient  of  species  z  in  layer  a 
X  number  of  water  molecules  per  sulfonic  acid  group 

in  Nation 

vpi  stoichiometric  coefficient  of  species  z  in  reaction  p 
vpe-  stoichiometric  coefficient  of  electrons  in  reaction  p 
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vPte-  stoichiometric  coefficient  of  electrons  in  the  RDS  in 

electrode  reaction  p 

§  electro-osmotic  drag  coefficient  of  water 

pM<A  anode  catalyst  metal  density,  g  cm-3 
pM<c  cathode  catalyst  metal  density,  g  cm-3 
crB  effective  PEM  conductivity 

X[  residence  time  of  species  i  in  electrode  compart¬ 

ment,  s 

(pi  fraction  of  metal  surface  in  contact  with  ionomer 

O  electrode  potential,  V 

OP)0  half-cell  thermodynamic  (reversible)  potential  of 
electrode  reaction  p,  V 
Oc  Cathode  potential,  V 

Oyi  o  Thermodynamic  potential  of  anode,  V 

oo Ru  mass  fraction  of  Ru  in  anode  catalyst 

Subscripts/superscripts 

rate-determining  step 
active  metal  catalyst  area  basis 
0  equilibrium 

A  anode  layer 

B  PEM  layer 

b  bulk-phase 

C  cathode  layer 

D  anode  diffusion  layer 

E  cathode  diffusion  layer 

e-  electron 

/  interfacial  or  ionomer 

i  species 

in  feed 

L  limiting 

Me  methanol 

ref  reference 

W  water 

X  crossover 

a  layer  in  MEA 

p  electrode  (anode  or  cathode) 


proved  challenging  because  of  multiple  nonlinearities.  Wang  and 
Wang  [17]  devised  a  comprehensive  two-phase  multicomponent 
computational  fluid  dynamics  (CFD)  model  capable  of  predicting 
the  mixed  potential  effect  at  the  cathode  from  methanol  crossover, 
as  well  as  the  complete  polarization  plot.  Chen  and  Zhao  [18]  devel¬ 
oped  a  model  for  a  passive  DMFC  to  investigate  the  effects  of  heat 
and  mass  transport  and  electrochemical  kinetics  on  overall  per¬ 
formance.  Rice  and  Faghri  [19]  investigated  the  one-dimensional 
mass  transport  of  methanol  vapor,  water  vapor,  and  carbon  diox¬ 
ide  in  a  vapor-fed  passive  DMFC  in  order  to  numerically  describe 
its  temporal  performance. 

On  the  other  hand,  there  are  simpler,  semi-empirical,  alge¬ 
braic  models  available  in  the  literature,  which  are  convenient 
for  real-time  simulation  and  control.  Thus,  Chiu  et  al.  [4,20] 
have  developed  an  algebraic  semi-empirical  model.  Using  the 
Levenberg-Marquardt  algorithm,  8  model  parameters  were  fitted 
based  on  regressing  a  set  of  28  measurements.  Using  an  additional 
1 12-measurement  data-set  for  model  verification,  they  found  that 
the  proposed  model  was  able  to  accurately  predict  methanol 
crossover  as  well  as  the  DMFC  performance  and  efficiency  over 
a  broad  range  of  conditions.  Dohle  and  Wippermann  [21]  devel¬ 
oped  a  semi-empirical  model  describing  the  relationship  between 
current  and  anode  overpotential,  cell  potential  versus  current,  and 
methanol  permeation. 

This  paper  provides  an  isothermal  model  of  the  DMFC  that  is 
based  on  a  detailed  consideration  of  the  key  transport  and  reaction 
processes,  but  uses  defensible  simplifying  assumptions  to  obtain 
analytical  solution.  Thus,  explicit  analytical  expressions  are  derived 
for  the  anode  and  cathode  overpotentials  in  terms  of  current  den¬ 
sity,  exchange  current  density,  the  limiting  current  density,  and  the 
crossover  limiting  current  densities,  in  turn  described  in  terms  of 
basic  mass  transfer  and  kinetic  parameters.  Further,  since  oxygen 
can  also,  in  principle,  diffuse  over  to  the  anode  in  a  DMFC  [22],  or 
can  infiltrate  from  the  CO2  vent  in  a  passive  DMFC,  the  potential 
effect  of  any  oxygen  present  at  the  cathode  is  also  included.  Addi¬ 
tionally,  power  and  heat  generated  in  a  DMFC  is  evaluated  along 
with  its  efficiency.  The  model,  involving  a  priori  parameters,  with 
the  exception  of  two  fitted  parameters,  is  able  to  accurately  predict 
the  comprehensive  set  of  experimental  data  of  Chiu  et  al.  [4,20] 
over  a  broad  range  of  conditions. 


concomitantly  reduces  methanol  crossover  [11].  Finally,  higher 
flow  rates  at  anode  and  cathode  can  significantly  reduce  mass 
transfer  limitations  [12]  and  more  effectively  remove  CO2  and  H20 
produced  at  the  anode  and  cathode,  respectively. 

Over  the  past  decade,  there  has  been  a  significant  advancement 
in  the  theoretical  understanding  of  these  issues.  Only  a  hand¬ 
ful  of  investigations  are  cited  here  for  brevity.  Thus,  Sundmacher 
and  Scott  [13]  developed  a  multilayer  model  for  the  membrane- 
electrode  assembly  (MEA),  including  unsteady-state  mass  balance, 
charge  balance,  diffusion,  and  electrode  reaction  kinetics,  as 
well  as  thermodynamics  of  the  two-phase  C02 -methanol-water 
gas-liquid  mixture  phase  equilibria  in  the  anode-diffusion  layer 
(ADL).  However,  they  neglected  mass  transfer  limitations  in  the 
catalyst  layer  (CL)  because  of  its  thinness  compared  to  the  other 
layers.  Garcia  et  al.  [14]  included  the  CL  diffusion  limitations  in 
their  semi-analytical  solution  approach.  Their  results  indicate  that 
the  diffusion  limitations  are,  in  fact,  generally  small  in  the  anode  CL. 
Assuming  Tafel  kinetics,  Kulikovsky  [15]  solved  the  ACL  diffusion 
problem  analytically,  but  only  for  asymptotically  large  and  small 
currents. 

Meyers  and  Newman  [16]  developed  a  detailed  theoretical 
framework  for  describing  non-ideal  multicomponent  transport 
in  the  various  layers  of  the  MEA,  along  with  a  multistep 
MOR,  and  used  a  finite-difference  scheme  to  solve  the  cou¬ 
pled  differential-algebraic  equation  system.  However,  the  solution 


2.  Theory 

2.1.  Model  description 

Fig.  1  schematically  shows  the  concentration  profiles  of 
methanol,  oxygen,  carbon  dioxide,  and  water,  as  well  as  the  anode 
and  cathode  potentials,  <&A  and  Oc,  across  the  various  layers  of  the 
DMFC  MEA. 

Due  to  mass  transfer  limitations  in  the  bulk  phase  of  anode 
chamber  as  well  as  in  the  ADL,  layer  D,  thus,  the  concentration 
of  methanol  drops  as  it  is  consumed  via  MOR  in  the  anode  cata¬ 
lyst  layer  (ACL),  layer  A,  to  produce  C02,  protons,  and  electrons. 
Furthermore,  methanol  concentration  drops  across  the  PEM,  layer 
B,  as  methanol  is  transported  across  it  via  diffusion  and  electro- 
osmotic  drag  to  arrive  at  the  cathode,  where  it  also  undergoes 
MOR.  The  protons  produced  at  the  anode  diffuse  through  the  PEM 
to  the  cathode,  whereas  electrons  commence  their  journey  across 
the  external  circuit,  arriving  at  the  cathode,  depleted  from  hav¬ 
ing  performed  external  electrical  work.  Unlike  as  shown  in  the 
schematic,  however,  the  methanol  concentration  actually  declines 
somewhat  within  the  ACL,  as  it  is  consumed  in  the  MOR,  and  should 
properly  be  analyzed  via  a  distributed  diffusion-reaction  model 
[14,15].  However,  in  order  to  obtain  a  simple  analytical  solution, 
it  is  assumed  here  to  be  uniform  [13],  based  on  considering  the  ACL 
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Fig.  1.  Schematic  representation  of  a  DMFC  MEA,  including  concentration  profiles  shown  here  for  a  concentrated  methanol  feed,  external  current  and  internal  crossover 
currents,  and  potential  profile  across  the  various  MEA  layers. 


as  a  thin  layer,  which  is  also  tantamount  to  assuming  a  small  Thiele 
modulus. 

The  MOR  reaction,  which  is  very  sluggish  owing  to  the  strong 
CO  catalyst  adsorption  at  low  temperatures,  is  facilitated  by  the 
anode  overpotential,  r]A  =  <&A-  ®A  0,  which  in  turn  is  determined 
by  the  drawn  anode  current  density,  iA,  comprising  of  both  the 
external  current  i,  as  well  as  any  current  iA>x,  at  the  anode  due  to 
crossover/presence  of  02  in  the  anode  chamber,  i.e.,  iA  =  i  +  iA<x.  The 
C02  that  is  produced  at  the  anode  can  either  diffuse/bubble  across 
the  ADL  to  be  vented,  or  can,  due  to  its  high  permeance  [23],  also 
diffuse  across  the  PEM  to  be  removed  partly  on  the  cathode  side. 
Of  course,  a  considerable  amount  of  C02  is  also  produced  at  the 
cathode  from  methanol  crossover. 

Similarly,  oxygen  in  the  cathode  chamber  diffuses  across  the 
cathode  diffusion  layer  (CDL),  layer  E,  to  arrive  at  the  cathode  cat¬ 
alyst  layer  (CCL),  layer  C,  where  it  undergoes  the  oxygen  reduction 
reaction  (ORR)  (Fig.  1).  The  ORR  involves  protons  and  electrons 
produced  at  the  ACL  and  arriving  via  PEM  and  the  external  cir¬ 
cuit,  respectively.  In  addition,  however,  protons  and  electrons 
produced  in  situ  at  the  CCL  via  MOR  are  utilized  in  the  ORR. 
Thus,  the  cathode  current  density,  ic  =  i  +  icj( •  The  cathode  over¬ 
potential,  ric=<&c-®c}o,  that  is  negative  (^7c < 0),  accelerates  the 
ORR,  which  is  otherwise  exceedingly  sluggish.  The  CCL  is  also 
considered  as  thin,  i.e.,  without  diffusion  limitations,  so  that 
the  concentrations  within  it  are  assumed  to  be  uniform  as  well 
(Fig.  1). 

The  methanol  arriving  at  the  cathode  also  undergoes  elec¬ 
trochemical  MOR,  rather  than  chemical  oxidation,  as  sometimes 
believed  [24],  due  to  the  huge  overpotential  there  for  the  MOR, 
*M,x  =  -  ^a.o  (Fig.  1),  even  when  the  external  current  i  =  0,  i.e., 

under  open  circuit  voltage  (OCV)  conditions,  effectively  creating 
an  internal  short  circuit,  or  a  crossover  current  (ic,x)-  In  fact,  since 
the  thermodynamic  voltage,  V0  =  Oc,o  -  ^a.o  [2],  we  may  alter¬ 
nately  write  i)A^  =  Vo  +  r]c,  as  shown  in  Fig.  1.  As  a  result,  we  will 
assume  below  that  the  concentration  of  methanol  in  the  CCL^  0. 
This  assumption,  as  well  as  small  concentration  variations  within 
the  CL,  is  supported  by  the  results  of  Garcia  et  al.  [14].  Likewise, 
any  oxygen,  albeit  to  a  much  smaller  extent,  permeating  from  the 
cathode  can  also  readily  undergo  electrochemical  ORR  due  to  large 
overpotential  r]C^  =  <&A  -  Oc,o,  thus  creating  an  internal  crossover 
current  (iAX). 


2.2.  Basic  constitutive  relations 


In  order  to  develop  a  simple  albeit  robust  analytical  model,  and 
as  is  not  uncommon,  the  electrode  kinetics  of  the  MOR  and  ORR, 
denoted  generically  as  the  electrode  reaction  p,  are  approximated 
here  by  the  Butler-Volmer  equation,  though  strictly  this  is  valid 
only  for  the  elementary  steps  involved  in  an  overall  electrode  reac¬ 
tion  [  2  ] ,  the  overall  electrode  kinetics  in  reality  being  more  nuanced 
[16,25] 


p 


7?  =  ~r  =  exP 
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{<Xp)Vpe-Fllp 
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RT 


■  j  -  exp  j- 


(oip)Vpe-Fllp  \ 

RT  J 
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where  i*  0  (or  ip0)  is  the  exchange  current  density,  vpe -  is  the  stoi¬ 
chiometric  coefficient  of  electrons  in  the  electrode  reaction,  and  ap 
and  bip  are  the  fitted  transfer  coefficients  of  the  forward  and  reverse 
electrode  reaction.  Here,  the  superscript  asterisk  denotes  current 
density  based  on  actual  active  electrocatalyst  surface  area  (ECSA), 
while  that  without  it  denotes  current  density  in  terms  of  geometric 
MEA  area. 

Alternately,  the  transfer  coefficients  may  be  construed  simply  as 
the  symmetry  factor  of  the  rate-determining  step  (RDS),  typically 
taken  as  1  \2  for  elementary  reactions.  In  such  an  event,  i.e.,  when 
it  is  justifiable  to  assume  that  ap=ap  =  ap,  where  the  superscript 
dot  implies  RDS,  the  Butler-Volmer  equation  can  be  cast  into  the 
convenient  form  [2] 
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RT 
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Here,  the  current  density  per  unit  MEA  area,  and  kinetics  of  the 
electrode  reaction  are  interrelated  via  [2] 


iP  =  Ym?p  =  YM(Fvpe-r*y,  iAo  =  yMi*pfi  =  VM(FVpe-r*0)  (3) 

as  are  the  exchange  current  density  and  the  corresponding  rate. 
Here  the  roughness  factor  (yM)  is  the  ratio  between  the  active 
electrocatalyst  area  (ECSA)  and  the  geometric  area  of  the  elec- 
trode/MEA.  Further,  the  rate,  rp,  defined  on  the  basis  of  the 
geometric  MEA  area,  rp  =  Ym^p- 

The  electrode  (both  anode  and  cathode)  reaction  rate  is  further 
assumed  to  be  pseudo  first-order  [2],  along  with  the  treatment  of 
the  electrode  layer  as  an  interface,  i.e.,  a  layer  of  small  thickness, 
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so  that  the  concentration  of  a  species  i  within  it  is  uniform  (Fig.  1), 
e.g.,  for  the  anode,  citA  =  ci<A{a)  =  ciiA(b ),  leading  to 


rp  =  k*pcUA{ay,  r*0  =  k*p<t,oKUAcUb  (4) 

where  ci  b  is  bulk-phase  concentration  of  species  i  in  the  electrode 
chamber  or  flow  channels,  and  Ki<A  is  the  partition  coefficient  of  the 
species  i  from  the  bulk-phase  to  the  ACL,  defined  below. 

More  correctly,  of  course,  r*  =  spkpCj(a),  in  Eq.  (4),  where  ep  is 
the  effectiveness  factor,  e.g.,  sp  =  tanh  4>p!4>p,  where  4>p  is  the  Thiele 

modulus  for  the  electrode  reaction,  4>p  =  L^Jk^/D^ ~,  and  L  is  the 

electrode  layer  thickness.  However,  numerical  results  [14]  indicate 
that,  for  the  relatively  low  current  densities  typically  encountered 
in  DFMCs,  the  assumption  of  unit  effectiveness  factor  is  acceptable, 
as  implicit  in  Eq.  (4)  above. 

For  the  one-dimensional  (z-direction)  flux  (moles  per  unit  area 
of  the  MEA)  of  a  species  z,  the  Nernst-Planck  equation  is  used  [26] 


N;,z  =  - 


De(^i+Z.C.L^) 

1  \dz  +  '  'RT  dz ) 


+  xtN 


(5) 


where  the  terms  in  the  parentheses  represent  diffusion  and  electro¬ 
migration,  respectively,  while  the  last  term  represents  convection. 
The  electromigration  term  is,  of  course,  zero  for  uncharged  species 
(z;  =  0).  Further,  the  relationship  between  flux  of  a  species  partici¬ 
pating  in  the  electrode  reaction  and  the  electrode  current  density 
is 


ip  =  Ni,z 


(6) 


where,  e.g.,  the  stoichiometric  coefficient  of  methanol  in  the  MOR, 

VA,  Me  =  -  1  • 

The  concentration  of  a  species  z  in  a  layer  a  of  the  MEA  in  equi¬ 
librium  with  the  bulk  phase  concentration  is  described  in  terms  of 
the  partition  coefficient 


^i,a 

Ci,b 


(7) 


which  is,  of  course,  a  function  of  temperature  and,  for  non-ideal 
mixtures,  of  concentration.  Pseudo-equilibrium  is  assumed  at  the 
interface  between  two  phases/layers.  For  instance,  the  partition 
coefficient  of  methanol  [27]  from  bulk  to  the  ADL,  /cMe,D»  is  deter¬ 
mined  in  part  by  the  fraction  of  the  pore  volume  occupied  by  gas 
bubbles,  and  the  vapor-liquid  equilibrium  [13].  In  the  absence  of 
a  gas  phase  within  the  ADL,  it  is  approximately  unity  for  a  liquid 
feed,  but  not  so  for  a  vapor  feed. 

Across  the  bulk  in  the  anode  chamber  as  well  as  in  the  adjacent 
gas-diffusion  layer  (GDL),  flux  of  a  species  z  is  constant  at  steady 
state,  i.e., 


D  —  K  \  ^i.b 


Q,d(0) 

Ki,D 


DId 

Ld 


{q,d(0)  -  Q  D(a)} 


(8) 


where  the  first  equation  represents  the  flux  from  the  bulk  phase 
(e.g.,  liquid  or  vapor  phase  methanol  feed  in  the  anode,  or  gas-phase 
02  in  the  cathode)  to  the  GDL,  while  the  second  equation  is  the  flux 
through  the  GDL.  These  two  can  be  combined  into  the  form 


NUD  =  PftDUub- 


cuda) 

kld 


O) 


A  similar  relation  may  be  written  for  the  cathode  diffusion  layer 
(CDL),  layer  E. 


2.3.  Methanol  crossover  and  anode  overpotential 


Methanol  that  escapes  reaction  within  the  ACL  can  permeate 
readily  through  the  PEM  (layer  B  in  Fig.  1 )  via  diffusion  and  electro- 
osmotic  drag 

^Me,B  =  |cMe,B(^)  -  CMe,B(c)  j-  +  —  —j-  ^  ^  CMe.fidZ  (H) 

represented,  respectively,  by  the  first  and  second  terms  on  the 
right-hand  side  of  this  equation.  This  equation  results  from  the 
integration  of  Eq.  (5)  for  constant  flux,  in  the  absence  of  electro¬ 
migration,  and  with  the  assumption  that  the  methanol  mole  frac¬ 
tion  within  the  PEM  layer,  XMe,B  =  CMe,Bl{CMe,B  +  Cw,B)^CMe,BlCw,B, 
because  of  the  high  molarity  of  water  as  compared  to  methanol. 
Here,  §  is  the  electro-osmotic  drag  coefficient,  roughly  2.5  -  3.5 
molecules  of  water  per  proton  transported  across  the  PEM  [28]. 

As  assumed  above,  the  concentration  within  the  ACL, 
cMeA  =  cMeAa)  =  cMeA(b),  and  further,  that  any  methanol  crossing 
over  is  consumed  entirely  within  the  CCL,  i.e.,  cMe,B(c)  ->  0.  These 
simplifications,  along  with  the  assumption  of  linear  variation  of 
methanol  concentration  within  the  electrolyte  layer,  B,  reduce  the 
membrane  flux,  Eq.  (11),  to 


NMe,B 


DM  e,B 

Lb  F  2  cw,b  J 


(  XMe,B 
\XM  e,D 


Qvie,a(a) 


(12) 


where  it  has  been  further  assumed  that  the  partition  coefficient  of 
methanol  in  the  ACL  is  the  same  as  that  in  the  ADL,  i.e.,  /cMe^  =  /cMe,D- 
Of  course,  using  Eq.  (6)  with  this,  the  crossover  current  due  to 
the  methanol  flux  across  the  PEM  is 


Me 


VA,e- 
~VA,  Me 


2  Cw,B 


/  ^Me,B 
\XM  e,D 


Cm  e,a(a) 


(13) 


We  will  find  it  useful  to  define  a  limiting  methanol  crossover  cur¬ 
rent  density,  i.e.,  the  largest  possible  methanol  diffusive  flux  across 
the  PEM,  which  is  under  OCV  conditions,  i.e.,  when  the  external 
current,  z  =  0.  Further,  we  imagine  that  under  these  conditions, 
there  is  little  diffusion  limitation  across  the  ADL  either,  so  that 
CMe^(a)  =  CMe,D(a)  =  cMe,D(0)  =  /CMe,DCMe,b  (Fig-  1).  Fe..  the  methanol 
bulk  concentration  is  available  at  the  ACL-PEM  interface.  Thus,  the 
corresponding  limiting  crossover  current  density  of  methanol,  from 
Eq.  (13),  is 


*X,M  e,L  = 


De 

uMe,B 

Lb 


VA,e~ 
~VA,  Me 


^Me,£?CMe,b 


(14) 


In  reality,  the  crossover  current  at  OCV  (i.e.,  when  z  =  0)  would  be 
lower  due  to  the  finite  methanol  mass  transfer  resistance  across 
the  ADL,  as  discussed  below.  At  any  rate,  using  this  definition,  Eq. 
(13)  is  rewritten  as 


lx,  Me 
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where  PeD  is  the  effective  (overall)  permeance,  or  effective  (over¬ 
all)  mass  transfer  coefficient,  including  both  bulk  mass  transfer 
coefficient  as  well  as  GDL  permeance,  Pl  D 


1  _  1  1  p  _KUDDeiD 

Pf,D~Pi,D  k,'  UD~  LD 


(10) 


To  determine  this  crossover  current,  thus,  cMe<A{a )  is  needed.  For 
this,  we  begin  with  a  mass  balance  written  around  the  ACL, 
assumed,  as  mentioned  above,  to  be  lumped,  along  with  the  use 
of  Eq.  (4),  to  provide 


NMe,D(a)  -  NMe,s(b)  =  (— V/\,Me)Ol  =  (-VA,Me)yMfe4CMe,/i(a)  (16) 
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This  equation  can  be  used  to  solve  for  cMe<A{a)  by  using  in  it  Eqs. 
(9)  and  (12),  along  with  CMe,D(a)  =  cMe^(a),  since  the  partition  coef¬ 
ficients  in  ADL  and  ACL  have  been  assumed  to  be  the  same,  i.e., 


CMe,A(a) 

K'Me^Qvie.b 

1  +  {(^Me^/^8)  +  O'/*7 )(£/2cW,b)}  (KMeV^Me,!))  +  {~vA,Me)YMk^KMe,D/P^ieD 

(17) 

where  the  last  term  in  the  denominator  represents  the  ratio 
of  the  maximum  rate  of  reaction  within  the  ACL  (when 
cMeA  =  cMeAa)  =  l<:Me,DCMe,b)  t°  the  maximum  rate  of  diffusion 
through  the  ADL  (when  cMe^(a)  -*  0),  while  the  second  term  in  the 
denominator  represents  methanol  crossover. 

Before  returning  to  the  crossover  current,  let  us  first  use  this 
expression  in  Eq.  (3),  along  with  Eq.  (4),  to  provide  the  anode  cur¬ 
rent  density 


.  YmFVaq-  k^K:Me,DcMe,b 

1  +  {(*W«  +  (i/f X?/2Cw,B)}  (KMe,B/PeMe,D)  + 

(18) 

When  the  rate  of  diffusion  in  the  ADL  is  significantly  less  than  the 
rate  of  the  overall  electrode  reaction,  the  “1”  may  be  neglected  com¬ 
pared  to  the  other  terms  in  the  denominator.  Of  course,  then  the 
crossover  term  is  also  negligible,  as  then  CMe,D(a)-^  0.  This,  there¬ 
fore,  defines  the  ADL  diffusion-restricted  anode  limiting  current 
density  from  above,  i.e., 

*A,L  =  ( — 77  ^  ^Me,DCMe,b  (19) 

y  —  PA, Me  J 

corresponding  to  the  maximum  possible  diffusion  flux  across  the 
ADL. 

The  previous  two  equations  can  be  combined  and  rewritten  as  a 
ratio  between  the  anode  current  density  and  its  diffusion-restricted 
limiting  current  density 
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(20) 

which  can  be  further  rearranged  as 
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Using  Eqs.  (14)  and  (19)  for  the  term  within  the  square  brackets  on 
the  right-hand  side  of  this  equation,  it  may  be  recast  as 
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where  the  ratio  of  concentrations,  cMetb/cWtB  =  (cMetb/cWtb) 
(cwf/cw^^XMe.b/Kwj,  where  xMe>b  is  the  mole  fraction  of 
methanol  in  the  bulk  and  /cw,b  is  the  partition  coefficient  of  water 
in  the  membrane,  i.e.,  /cw,b  -  Cw,B/cw,b 

Also,  from  Eqs.  (3),  (4)  and  (19),  we  can  write  ratio  of  the  anode 
exchange-current  density  to  its  diffusion-restricted  limiting  cur¬ 
rent  density  as 


iA  o  (-vA,Me)YM^A  0QKMe,D 
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Combining  the  last  two  equations,  thus,  the  ratio  of  anode  rate 
constant  to  that  at  equilibrium  anode  potential 
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Finally,  combining  this  with  Eq.  (2)  provides 
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which  may  be  further  inverted  to  provide  an  explicit  equation  for 
the  anode  overpotential 


0a  = 


RT 


UAVA,e-F 


-sinh 
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being  the  combined  overpotential  associated  with  kinetic  and  dif- 
fusional  processes,  both  in  the  bulk  and  in  the  GDL,  as  well  as 
the  effect  of  methanol  crossover  on  it.  In  the  above,  we  have  used 
VA,e-/{-VA, Me)  =  +6/{-(-l)}  =  +6  for  the  MOR. 

Returning  now  to  the  issue  of  methanol  crossover  current,  using 
Eqs.  (14),  (17)  and  (22)  in  Eq.  (15), 

{*X,Me,L  +  3(/CMe,B/^W,B)^Me,b?(0}  (1  -  U^A,l) 

*X,Me  —  - ; - y~. - (27) 

[1  +  (1/U,l)  {*X,Me,L  +  3(/CMe,B/^W,B)^Me,b?(0}] 

where  again  v^?e-/(-^,Me)  =  6  for  MOR  has  been  used,  and  we 
have  further  assumed  in  the  last  term  in  the  numerator  that  the 
anode  current,  iA  =  i  +  ix,o2  ~  i,  i.e.,  ix,o2  ->  0. 

Finally,  under  OCV  conditions  (i  =  0),  this  reduces  down  to 
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(28) 


which,  as  alluded  to  above,  suitably  tempers  the  limiting  crossover 
current  density  defined  by  Eq.  (14)  by  the  diffusion  limitations 
within  the  ADL  under  OCV  conditions. 


2.4.  DMFC  polarization 


The  observed  voltage  in  a  DMFC  is  determined  by  the  difference 
between  the  thermodynamic  voltage  and  the  anodic  overpotential, 
the  cathodic  overpotential,  the  Ohmic  losses  in  the  PEM,  as  well  as 
any  internal  resistance  exhibited  by  the  PEM  [2,22],  i.e., 

v  =  v0  -  rjA  +  ric  +  riB  -  Vi  (29) 


where  the  anode  polarization,  with  the  anode  current,  iA  =  i  +  tx,o2  * 
in  Eq.  (26)  is 
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i.e.,  the  anode  overpotential  is  affected  by  both  methanol  crossover 
as  well  as  oxygen  crossover  or  alternate  infiltration  into  the  anode 
chamber.  This  effect  of  the  presence  of  oxygen  at  the  anode  on  its 
polarization  has  not  so  far  been  discussed  in  the  literature.  Neither 
has  the  effect  of  methanol  crossover  on  anode  overpotential. 
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In  a  manner  similar  to  that  described  in  the  last  sub-section  for 
anode,  the  cathode  polarization  may  also  be  determined,  providing 


Finally,  the  power  density  (W  cm-2  MEA)  produced  by  the  DMFC 
can  be  obtained  simply  from  P= iV,  i.e., 
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where  it  is  assumed  that  the  effect  of  the  crossover  of  oxygen 
through  the  membrane  on  cathode  polarization  is  small,  i.e.,  it  can 
be  treated  virtually  as  impervious  [2,22],  so  that  the  cathode  over¬ 
potential  is  affected  by  only  methanol  crossover.  Further,  r/c  <0 
here,  of  course,  because  ”ce-  <  0  [2], 

The  CDL  diffusion  restricted  cathode  limiting  current  density,  in 
analogy  to  Eq.  (19)  is 


ic,L^4FPe02Ec0^b  (32) 

while  the  effective  permeance  of  oxygen  in  the  CDL  is  given  by  Eq. 
(10),  i.e., 
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The  oxygen  crossover  current  may  be  roughly  estimated  as  follows. 
There  are,  in  fact,  two  sources  of  oxygen:  (1)  in  anode  methanol 
feed  itself,  assuming  it  is  saturated  with  oxygen  from  the  air,  and 
(2)  from  the  cathode  air  feed.  Assuming  that  the  anode  oxygen  at 
the  ACL  is  limited  by  the  diffusion  across  the  ADL,  Eq.  (19)  can  be 
adopted  for  it,  while  the  cathode  oxygen  is  limited  by  diffusion 
across  the  CDL  as  well  as  the  PEM  layer,  which  can  be  approximated 
using  Eq.  (28).  Thus,  assuming 


ix,02  =  |po2,D+(^  +  ^)  J  4FCq2  b 

^  (^02,D +^02,B)4^c02,b  (34) 

where  the  simplification  indicated  in  the  second  equality  in  Eq.  (32) 
results  because  the  permeance  of  02  in  CDL  is  much  larger  than  the 
permeance  of  O2  in  the  PEM.  This  is  because  while  the  diffusion  in 
the  CDL  is  mostly  in  the  gas  phase,  that  in  the  ADL  is  mainly  in  the 
liquid  phase. 

In  the  above,  the  oxygen  permeance  in  PEM 
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while  the  effective  permeance  of  oxygen  in  the  ADL  is  given  by  Eq. 

(10). 

Finally,  the  combination  of  Eqs.  (29),  (30)  and  (31)  provides  the 
DMFC  polarization 


1/  1/  ■  u-i 

V  =  V0-  — - -sinh 

aAvA,e-F 


{i  +  ix,o2  )/to, 
!-(*'  +  ix,  o2 


to,  0  ) 

)/to,r  J 

{1  +  i(Wl+3(S)Wi)}_ 

sinh-1  [U  (1  +  ,^e)/lC0  \  1 
F  \2  \l-(*  +  *X,M e)/lC,L  J  J 


acvc,e- 


iLb 

&B 


-  KRi ) 


(36) 


where  the  second-to-last  term  on  the  right-hand  side  is  the  Ohmic 
overpotential  in  the  electrolyte  layer,  and  1?/  is  any  interfacial  MEA 
resistance  due  to  imperfect  bonding  of  different  layers  or  due  to 
partial  delamination  with  use. 


P  =  i Vq  -  .  iRT  sinh 
aAvA,e-F 


1  if  (i  +  ix,  c 
2  \  1  -(i  +  l 


o2)/to,o 


ix,  o2  )/to,  j 


K(‘ 


x  <  1  +  1 -  (  ix,Me,L  +  3  [  'Me,B  ]  XMe 


iRT  .  u. 

+  - — 7 - -sinh 

acVC,e-F 

®B 


,e.)} 

-1  f  1  f  (*  +  *X,Me)/*C,0  \  1 

\  2  \  1  -  (i  +  ix,Me)/ic,L  J  J 


(37) 


This  completes  the  development  of  our  analytical  model  for  the 
DMFC.  It  is  seen  that  it  is  quite  tractable  and  that  many  of  the  terms 
in  the  equations  have  tangible  meaning. 


2.5.  Model  parameters 


The  veracity  of  above-developed  analytical  model  and  its  under¬ 
lying  assumptions  is  ascertained  by  comparing  it  to  experimental 
data,  using  the  best  available  model  parameters  in  the  literature, 
with  a  minimal  of  fitted  parameters. 

The  thermodynamic  DMFC  cell  potential  in  Eq.  (36)  is  given  by 


RT 

V0  =  1.214-1.4  x  10“4(r-298)+  — In 

or 


aMe 

<3H20 


(38) 


where  the  first  term  on  the  right-hand  side  represents  the  stan¬ 
dard  DMFC  cell  voltage,  V°  =  1.214  V,  corrected  for  temperature 
and  species  activities  by  the  remaining  two,  with  partial  pressures 
given  in  bars  or  in  atmospheres.  We  assume  the  activity  of  water  to 
be  unity,  and  replace  the  activity  of  methanol  by  its  mole  fraction 
for  dilute  feeds.  Then,  under  the  usual  DMFC  operating  conditions, 
the  thermodynamic  DMFC  cell  potential  V0  « 1.18  V. 

Let  us  next  consider  electrode  kinetics.  Both  the  MOR  [29]  and 
the  ORR  [30]  are  very  complex  electrode  reactions  with  multiple 
steps  and  pathways.  There  is  no  definitive  word  yet  in  the  literature 
on  the  mechanism  or  kinetics  of  either  of  these  important  reactions. 
In  fact,  even  an  unequivocal  value  for  the  exchange  current  density 
of  MOR  is  not  yet  available. 

For  the  ORR,  as  usual,  we  assume  the  Butler-Volmer  form  of  Eq. 
(2)  along  with  ac  =  1/2  and  vc,e-  =  -2,  construed  as  those  corre¬ 
sponding  to  the  RDS,  where  the  exchange  current  density  is  given 
by  [22] 


ic,  0  =  Ym,c 


C,0,ref 


(39) 


where  c02,o  is  the  concentration  in  the  CCL  under  equilibrium 
conditions,  i.e.,  c02,o  =  Ko2,cCo2,b>  als°  accounting  for  the  water 
vapor  pressure,  along  with  the  reference  exchange  current  density 
*'c  0  ref  =  1-0  x  10-10  A  cm-2  at  Tref=  298  K,  reference  concentra¬ 
tion  corresponding  to  Po2,re/  =  1.0  atm,  and  an  effective  activation 
energy  Ec  a>0  =  67,  000 J  mol-1  [22]. 

Here,  the  relationship  between  the  actual  rates  of  electrodic 
reactions  on  the  basis  of  unit  triple-phase  interface  (TPI)  area,  and 
those  observed  experimentally  on  the  basis  of  unit  geometric  MEA 
area  is  governed  by  the  electrocatalyst  surface  roughness,  which 
may  be  estimated  from  [2] 


(40) 


where  pm  is  the  catalyst  density,  771m  is  the  catalyst  leading,  dM 
is  the  supported  or  unsupported  catalyst  crystallite  diameter,  and 
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(pi  is  the  fraction  of  metal  catalyst  surface  in  contact  with  the 
ionomer,  forming  the  TPI.  It  may  be  noted  that  for  Pt-Ru  alloy, 
Pma  =  (1  -<^Ru)ppt  +  <^RuPRu  where  &>ru  is  the  mass  fraction  of  Ru, 
Ppt  =  2 1 .4  g  cm-3 ,  pRu  =  1 2 .3  g  orr3 . 

For  the  case  of  the  MOR  on  the  conventional  Pt-Ru  anode, 
the  dual  site  mechanism  is  broadly  accepted  [31],  wherein  the 
methanol  adsorbs  and  dissociates  into  adsorbed  carbon  monoxide, 
CO  •  Pt,  on  Pt  sites  [32],  covering  these  almost  completely  at  lower 
temperatures  and  overpotentials  because  of  its  strong  adsorption, 
while  water  dissociates  more  readily  on  the  Ru  sites  into  adsorbed 
OH,  i.e.,  OH  •  Ru,  which  finally  oxidizes  the  adsorbed  CO  into  C02, 
releasing  the  Pt  sites  on  which  the  electrochemical  MOR  can  then 
be  repeated,  starting  with  the  adsorption  of  methanol,  and  then 
proceeding  through  step-wise  dehydrogenation  to  produce  CO  •  Pt. 

Assuming  then  the  first  step  in  the  MOR  sequence, 
CH3OH  +  Pt  5=*  CH3OH  •  Pt,  to  be  the  RDS  in  the  sequence  [5],  we  may 
write  the  rate  of  the  electrode  reaction  as  ?!  «  /<iCMe(l  -0co-s)» 
in  which  6C o  s  represents  the  catalyst  site  fraction  occupied  by 
carbon  monoxide,  the  site  fractions  of  other  adsorbed  species 
(H  S  and  CH3OH  S,  etc.)  being  considered  negligible.  Thus,  the 
Butler-Volmer  form  of  Eq.  (2),  along  with  aA  =  1  /2,  and  i>Ae-  =  +1 
[33],  is  assumed  to  describe  the  potential  dependence,  where  the 
exchange  current  density  is  given  by 


to,  0  =  Vm,a 


cMe,0 
c  Me,  ref 


f  1  ~  0CO-S  A 

\  1  “  0 CO  S,  ref  J 


x  exp 


A,0,ref 


(41) 


where  cMe,o  is  the  ACL  concentration  under  equilibrium  condi¬ 
tions,  i.e.,  cMe,o  ^Me^Qvie^,  the  reference  methanol  concentration, 
cMe,re/=  ^-0  x  10-3  mol  cm-3,  and  the  MOR  activation  energy, 
EA  0q  =  65,  000J  mol-1  [34],  reported  for  Pt-Ru  catalyst  particles 
in  contact  with  Nation  electrolyte  [34].  This  is  close  to  the  value 
of  60,000Jmol_1  determined  by  Gasteiger  et  al.  [35]  for  liquid 
electrolyte,  and  is  further  consistent  with  first-principles  calcula¬ 
tions  [36].  In  comparison,  Gojkovic  et  al.  [33]  provide  an  activation 
energy  of  about  70  kj  mol-1 . 

The  fraction  of  sites  covered  by  CO,  $co  s>  iR  th e  exchange- 
current  density  expression  is  assumed  to  be  given  by  the  Langmuir 
isotherm 


Ocos  = 


^MeCMe,0 
1  +  ^MecMe,0 


(42) 


where  the  adsorption  equilibrium  constant  is  assumed  to  be  that 
for  CO  adsorption  on  Pt  [32] 


KMe  =  1.41  x  10  8  exp 


( 130,  000 J  mol 
^  RT 


(43) 


For  the  Pt-Ru  anode,  we  have  assumed  a  reference  intrinsic 
exchange  current  density  for  MOR  as  i*  0  ref  =  2  x  10-8  A  cm-2. 
This  is  one  of  the  two  main  fitted  parameters  in  our  model.  Although 
exchange  current  densities  for  MOR  have  been  reported  in  the  liter¬ 
ature  for  anode,  they  often  include  the  catalyst  surface  roughness, 
not  usually  reported.  In  other  words,  a  definitive  value  for  z*  0  re ^ 
has  not  so  far  been  reported. 

It  is  evident  from  this  that  even  on  state-of-the  art  Pt-Ru  cata¬ 
lyst,  the  MOR  is  an  exceedingly  slow  reaction,  much  as  ORR  is,  and 
is  a  key  reason  for  the  rather  lackluster  performance  of  DMFCs,  and 
a  requirement  for  high  catalyst  loading.  This  further  means  that, 
because  of  the  low  methanol  conversion  in  the  ACL  despite  the  high 
catalyst  loadings,  much  of  it  is  left  over  for  diffusion  across  PEM  to 
the  cathode,  further  impacting  the  cathode  polarization  and  also 
reducing  the  current  efficiency. 


Next,  let  us  consider  the  various  transport  parameters.  The 
methanol-water  binary  diffusion  coefficient  at  infinite  dilution  is 
fitted  to  the  available  experimental  data  [37]  as 

DMe,w  =  2.1  X  10-5  exp  j-20’/60  (i  -  3I3))  cm2  s_1  (44) 

The  effective  diffusion  coefficient  of  species  i  in  a  layer  a,  Deia ,  is 
dependent  on  the  volume  fraction  s  of  the  diffusing  phase  in  the 
layer  a  [2,13],  i.e.,  Deia  =  £^-5Di  a,  where  Dio(  is  the  molecular  diffu¬ 
sion  coefficient  in  phase  a ,  available  from  experiments  or  predicted 
from  theory  (as  in  gaseous  mixture),  and  the  exponent  1 .5  is  the  so- 
called  Bruggeman  exponent  of  volume  fraction  s ,  accounting  for 
both  the  porosity  and  the  tortuosity  factor  [38]. 

To  simulate  the  effect  of  carbon  dioxide  and  water  buildup 
within  the  ADL  and  CDL,  respectively,  on  the  effective  diffusion 
coefficients,  thus,  corresponding  volume  fraction  terms  ought  to 
be  included.  For  instance,  since  the  transport  of  oxygen  is  mainly 
in  the  gas-phase  due  to  its  low  solubility  in  water,  the  effective 
diffusion  coefficient  of  oxygen  within  the  CDL  is  [2] 

Deo2,E=45Cl-<lw)h5Do2,E  (45) 

where  the  effective  diffusion  coefficient  of  oxygen  Do2,e  in  the  mul¬ 
ticomponent  (N2,  H20,  C02)  gas-mixture  may  be  evaluated  from 
binary  diffusion  coefficients  as  described  in  [2],  and  qw  is  the  vol¬ 
ume  fraction  of  water  in  the  gas  pore  space  in  the  CDL,  assumed  as 
0.5  here  [2]. 

Similarly,  the  effective  diffusion  coefficient  of  methanol  in  the 
anode  gas-diffusion  layer  (layer  D)  may  be  estimated  from  the 
relation  DeMe  D  -  <Jco2)1'5£>Me,w  I13!-  where  <?co2  repre- 

sents  the  volume  fraction  of  the  ADL  pores  that  are  occupied 
by  C02  bubbles,  which  would  be  a  function  of  current  density 
as  well  as  methanol  feed  flow  rate  [39].  However,  unlike  the 
case  of  oxygen  transport  in  the  CDL,  the  transport  of  methanol 
in  the  ADL  occurs  to  significant  level  in  the  gas-phase  as  well 
[17]  because  of  its  high  concentration  as  well  diffusivity  in  the 
gas  phase.  Thus  a  better  estimate  of  methanol  effective  diffusion 
coefficient  in  ADL  would  include  its  transport  in  both  phases, 
i.e.,  D^je  D  =  sjj5  {(1  -  qc02)1'5^)Me,W  +  (to2)1'5/(:Me,C02^)Me,C02 
where  /cMe,c o2  is  the  liquid/gas  partition  coefficient,  which  varies 
with  temperature  as  well  as  with  methanol  concentration  [13],  and 
f^Me,co2  is  the  gas-phase  diffusion  coefficient  of  methanol.  Addi¬ 
tionally,  since  the  variation  of  qc o2  with  current  density  i  and  with 
temperature  is  not  well  characterized,  for  simplicity,  we  use  instead 
the  following  empirical  expression  for  the  effective  diffusivity,  fit¬ 
ted  to  the  limiting  current  density  data  of  Chiu  et  al.  [20] 

DeMe,D  =  9-75  x  10-1  exp  (-3(Lp5)  (46) 

which  is  the  second  main  fitted  parameter  in  our  model. 

The  effective  diffusion  coefficient  of  methanol  in  the  PEM  is 
assumed  to  be  similarly  given  by 

^Me,B  =  5^Me,W  (47) 

where  the  methanol-water  diffusion  coefficient,  DMe,w»  is  given 
in  Eq.  (44).  This  provides  an  activation  energy  for  the  effective 
methanol  diffusion  in  the  PEM  that  is  consistent  with  the  exper¬ 
iments  of  Ren  et  al.  [40]. 

In  this,  the  volume  fraction  of  water  in  the  PEM  sB  depends  upon 
A,  the  number  of  water  molecules  sorbed  by  the  electrolyte  layer 
per  acid  site,  i.e.  [2], 

e*-XT7  (48) 

where  r  =  Vr/Vw  =  537/18  «  30,  is  the  ratio  of  membrane  partial 
molar  volume  of  polymer  electrolyte  to  that  of  water. 
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Finally,  the  effective  oxygen  diffusion  coefficient  in  PEM  in  Eq. 
(33)  is  estimated  from 

do2,b  =  45do2,w  (49) 

where  Do2,w  is  the  molecular  diffusion  coefficient  of  oxygen  in 
liquid  water. 

There  are  few  good  estimates  available  in  the  literature  of  the 
mass  transfer  coefficient  Iq  in  the  bulk  (Eq.  (8))  for  the  electrode 
chambers  in  a  DMFC.  The  situation  is  complex  because  of  two-phase 
flow  in  the  flow  channels  of  both  the  electrodes.  For  instance  at  the 
anode,  the  mass  transfer  coefficient  depends  upon  the  hydrody¬ 
namics  in  the  channel  as  well  as  the  generation  of  CO2  bubbles  at 
the  GDL  surface,  the  extent  of  which  depends  upon  the  current  den¬ 
sity,  i  [13].  The  enhanced  convection  by  the  bubbles  would  increase 
/<Me-  On  the  other  hand,  these  bubbles  also  reduce  the  effective  area 
for  liquid-phase  mass  transfer  from  the  bulk  to  the  GDL.  In  fact,  Xu 
et  al.  [12]  found  little  effect  of  current  density  on  the  mass  transfer 
coefficient. 

In  the  absence  of  better  information,  for  now  we  simply  adopt 
the  following  correlation  for  bulk-phase  mass  transfer  coefficient 
under  laminar  flow  in  a  square  cross-section  channeled  monoliths 
of  side  d  and  length  L  [41  ] 

l<  d  /  d\0,81 

Sh  =  =  2.696 ^1  +  0.139  ReSc-J  (50) 

Such  relations  are  sometimes  written  in  terms  of  Peclet  number, 
Pe  =  Re  Sc  =  vd/Df.  Further,  we  can  write  the  velocity  in  terms  of 
the  volumetric  flow  rate  in  electrode  chamber,  v  =  Vp /(rid2 ),  where 
n  is  the  number  of  parallel  channels  in  the  flow  field,  e.g.,  n  =  1 
for  a  single  serpentine  channel.  Thus,  the  mass  transfer  coefficient 
becomes 


ki  =  2.696^ 


1+0.139 


d2/r A 


0.81 


Dj 


(51) 


where  rz  is  the  mean  residence  time  of  species  i  in  the  electrode 
chamber,  and  the  volume  of  the  electrode  chamber  =  nd2L,  i.e.,  rz  = 
Vp/nd2L.  Consequently,  the  second  term  in  the  parentheses  may 
be  construed  as  a  ratio  of  hydrodynamic  dispersion  to  molecular 
diffusion.  This  relation  could  be  assumed  to  provide  approximate 
mass  transfer  coefficient  values  for  both  the  methanol  feed  (vapor 
or  liquid)  as  well  as  for  the  cathode  air  feed. 

For  a  passive  DMFC,  thus,  since  00,  the  mass  transfer  coeffi¬ 
cient,  Iq  =  2.696(Dj /d).  Further,  the  dimension  d  over  which  diffusion 
occurs  under  quiescent  conditions  can  be  large  in  a  passive  DMFC, 
so  that  Iq  could  be  quite  small  and,  hence,  bulk  phase  diffusion  resis¬ 
tance  could  be  significant.  In  fact,  the  key  difference  between  active 
and  passive  DMFCs  is  the  much  higher  mass  transfer  coefficient  for 
methanol  as  well  as  for  oxygen  in  active  DMFCs. 

At  slow  flow  rates,  furthermore,  the  average  concentrations  of 
methanol  and  oxygen  in  the  bulk  in  electrode  chambers  are  lower 
than  the  feed  concentrations  due  to  consumption  at  the  electrodes. 
In  the  above,  in  fact,  we  have  assumed  that  the  bulk  concentrations 
of  methanol  and  oxygen  in  the  respective  electrode  chambers  are 
known.  Flowever,  what  is  known  are  the  inlet  feed  concentrations. 
The  actual  average  bulk  concentrations  in  the  electrode  chambers 
may  for  simplicity  be  taken  as  the  mean  of  the  inlet  and  outlet 
concentrations  [12],  i.e., 


Cib  = 


Q,  in  +  Q,out 


(52) 


where  the  inlet  and  outlet  concentrations  are  interrelated  via  a 
mass  balance  on  the  electrode  chamber  considered  as  a  well-mixed, 
i.e., 


in  ci,out)  —  ANiz 


(53) 


where  A  is  the  MEA  geometric  area.  Combining  the  above  two 
equations  and  relating  the  species  fluxes  to  corresponding  current 
densities  (Eq.  (6))  provides  for  the  anode  [12] 


^  A(i  +  ix, Me) 

HVIe,b  —  cMe,in  • 

l2tVA 


(54) 


while  for  the  cathode 


C(D 2,b  ~  ^02,in 


Ax 

8  FVC 


(55) 


where  the  O2  concentration  must  be  corrected  for  the  presence  of 
water  vapor  as  well  (Table  1 ). 

Finally,  the  Ohmic  losses  through  the  PEM  are,  of  course,  pro¬ 
portional  to  its  thickness  and  membrane  conductivity  (Eq.  (36)). 
Although  more  complex  forms  are  available  in  literature  for  mem¬ 
brane  conductivity  [42],  we  assume  the  simpler  expression  of 
Thampan  et  al.  [38] 


(56) 


where  the  parameter  8  varies  between  0.6  and  5.5  [2,38].  In  the 
following  8  =  3.5  was  chosen  as  it  agrees  well  with  the  conductivity 
data  of  Ren  et  al.  [40]  in  DMFC. 

The  degree  of  dissociation  of  acid  sites,  or  in  the  above,  is 


„  (X  +  l)-  -V++  l)2  -4A(1  -  l/K.) 

<>= - 2(1 -1/If.) -  1571 

where  X  is  the  number  of  water  molecules  sorbed  by  the  electrolyte 
layer  per  acid  site,  and  the  dissociation  equilibrium  constant  for 
acid  sites  is 


I<a  =  6.2  exp 


-52,300  n  _  1  \ 
VT  “  298 / 


R 


(58) 


In  Eq.  (56),  further,  sB,o  is  the  percolation  threshold,  roughly  cor¬ 
responding  to  A0  =  1.8.  Assuming  a  liquid  feed,  we  assumed  the 
membrane  to  be  saturated,  so  that  X  =  20  [40]. 

The  model  parameters  discussed  above  are  summarized  in 
Table  1. 


3.  Results  and  discussion 

To  check  if  our  model  as  developed  above  with  mostly  a  priori 
parameters  is  a  good  representation  of  the  DMFC  performance,  its 
predictions  are  compared  next  with  literature  data.  A  stumbling 
block  in  such  a  comparison  is,  of  course,  finding  a  comprehen¬ 
sive  and  consistent  set  of  data,  because  of  the  wide  variability  of 
DMFC  performance  reported  in  the  literature.  Chiu  et  al.  [4,20]  have 
recently  reported  such  a  set  of  data,  which  they  kindly  provided  to 
us  in  tabular  form,  and  is  used  below  in  our  comparison.  The  gen¬ 
eral  experimental  conditions  of  Chiu  et  al.  [4,20]  are  summarized 
in  Table  2. 

3.1.  DMFC  polarization  and  power  density 

Two  of  the  variables  that  most  profoundly  affect  DMFC  per¬ 
formance  are:  (1)  the  feed  methanol  concentration,  and  (2)  the 
operating  temperature. 

The  effect  of  methanol  feed  concentration  on  the  DMFC  perfor¬ 
mance  is  shown  in  Fig.  2  along  with  the  data  of  Chiu  et  al.  [20], 
for  fuel  concentrations  varying  from  0.75  M  to  1.5  M  and  air  at  the 
cathode,  and  at  a  constant  operating  temperature  of  313  K.  Fig.  2(a) 
shows  the  polarization  plot  (Eq.  (36)),  while  the  corresponding 
power  density  (Eq.  (37))  versus  current  density  is  shown  in  Fig.  2(b). 
It  can  be  seen  that  the  model,  with  the  set  of  parameters  listed  in 
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Table  1 

DMFC  model  parameters. 


Parameter 

Value 

Units 

Reference/comment 

AG° 

-702.7 

kj  K_1  mol-1 

Standard  Gibbs  free  energy  change 

A  H° 

-726.8 

kj  K_1  mol-1 

Standard  enthalpy  change 

V° 

0 

-AG°/6F=  1.214 

V 

Standard  thermodynamic  potential 

V° 

v  max 

-AH°/6F=  1.255 

V 

Standard  maximum  potential 

PPt 

21.45 

gem-3 

Platinum  density 

PRu 

12.3 

gem-3 

Ruthenium  density 

dwiA .  djw.c 

2.7 

nm 

Typical  for  C  supported  Pt  [33] 

< Pi 

0.75 

- 

[2] 

lA,0,ref 

2.0  x  10“8 

A  cm-2 

Fitted  (the  intrinsic  exchange  current  not  available  in  literature) 

lC,0,ref 

o 

O 

X 

o 

A  cm-2 

[22] 

Qvie.re/ 

1.0  x  10“3 

mol  cm-3 

Standard 

Cw 

55.5  x  10“3 

mol  cm-3 

Dilute  methanol  solutions 

Tref 

298 

K 

Standard 

Po2,ref 

1 

atm 

Standard 

aA 

1/2 

- 

Typical  for  elementary  RDS 

VA,e~ 

+1 

- 

Agrees  with  Tafel  slope  of  Gojkovic  et  al.  [33] 

ac 

1/2 

- 

Typical  for  elementary  RDS 

^C,e- 

-2 

K 

[22] 

Va,e- 

+6 

- 

Stoichiometric  coefficient  of  electrons  in  MOR 

Pd.Me 

-1 

- 

Stoichiometric  coefficient  of  methanol  in  MOR 

£a,cd0 

65 

kj  mol-1 

[34] 

£c<t>o 

67 

kj  mol-1 

[2] 

I<Me 

1.5  x  10~8exp  (13°'°°0J  m°rt,1 

1 

Estimated  as  CO  adsorption  equilibrium  constant  on  Pt 

Lb 

178  xl0“4 

cm 

Nation®  117  thickness 

Ld 

260 x 10“4 

cm 

[22] 

Le 

260  xlO-4 

cm 

[22] 

£d,  £e 

0.65 

- 

Typical  ADL/CDL  porosity 

ko2 

0.144 

- 

[2] 

Km  e,D 

1.0 

- 

Assumed 

Km  e,B 

0.4 

- 

[40] 

£ 

2.9 

- 

[28] 

E* 

14 

kj  mol-1 

Water  viscosity  activation  energy 

X 

20 

- 

[40] 

8 

3.5 

- 

Fitted  to  conductivity  data  of  Ren  et  al.  [40] 

Qw 

0.5 

- 

Volume  fraction  of  water  in  CDL  pore  s 

f^Me.W 

2.1  x  10-5  exp  | _  2M60  _ 

3Tb)} 

cm2  s_1 

[37] 

De 

M  e,D 

9.75  x  lO”1  exp  (-^^) 

cm2  s_1 

Fitted  diffusion  coefficient  in  ADL 

D  o2,w 

7.4  xW~8T(40-68t 

/x(25.60,6) 

cm2  s_1 

Oxygen  diffusion  coefficient  in  liquid  water 

Do  2,e 

0.357  x  (^)  ' 

cm2  s_1 

Gas-phase  oxygen  diffusion  coefficient 

Pw 

exp  (l  1.676  -HIS) 

atm 

Water  vapor  pressure  in  cathode 

r 

Vb/Vw  =  537/18 

- 

[2] 

Ri 

0 

£2  cm2 

MEA  interfacial  resistance 

Table  1 ,  fits  the  data  very  well  throughout  the  voltage-current  den¬ 
sity  range  investigated.  In  other  words,  the  model  appears  to  be 
robust. 

Although  not  readily  discerned  because  of  the  scale  in  Fig.  2(a), 
the  OCV  as  well  as  the  cell  potential  V  at  low  current  densities  i  is 


actually  higher  for  lower  methanol  feed  concentrations,  because 
of  the  lower  polarization  of  the  anode,  as  well  as  lower  methanol 
crossover  and  the  resultant  lower  polarization  of  the  cathode.  How¬ 
ever,  at  higher  current  densities  (e.g.,  >50mAcm-2  in  Fig.  2),  the 
trend  reverses,  with  higher  cell  potentials  resulting  at  a  given 


(a) 


V(V) 


0.25 


i  (A/cm2)  f  (A/cm2) 

Fig.  2.  Comparison  of  theoretical  predictions  with  experiments  of  Chiu  et  al.  [20]  for  (a)  polarization  plots  and  (b)  power  density  versus  current  density  of  a  DMFC  operating 
at  313 1<  with  air  at  the  cathode  and  different  methanol  feed  concentrations  in  the  range  of  0.75  M  to  1.5  M.  Table  2  provides  the  other  experimental  conditions. 
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Fig.  3.  Comparison  of  theoretical  predictions  with  experiments  of  Chiu  et  al.  [20]  for  (a)  polarization  plots  and  (b)  power  density  versus  current  density  of  a  DMFC  operating 
with  1  M  methanol  feed  and  air  at  the  cathode  at  different  temperatures  from  303  to  333  K.  Table  2  provides  the  other  experimental  conditions. 


current  density  for  the  higher  methanol  feed 
concentrations.  This  is  by  virtue  of  higher  diffusion  limita¬ 
tions  in  the  ADL,  i.e.,  lower  limiting  current  densities,  iAL,  at  the 
lower  methanol  feed  concentrations.  Likewise  in  Fig.  2(b),  the 
curves  are  bunched  together  at  lower  current  densities  but  then 
diverge  at  higher  current  densities,  when  the  higher  methanol 
feed  concentrations  can  produce  more  power.  The  power  density 
then  peaks  before  declining  as  the  DMFC  performance  becomes 
dominated  by  diffusion  limitations  in  the  GDL  at  higher  current 
densities.  Thus,  there  is,  an  optimum  power  density  for  a  given 
methanol  feed  concentration. 

The  effect  of  temperature  on  DMFC  performance  is  illustrated  in 
Fig.  3,  where  polarization  plots  are  provided  in  Fig.  3(a)  at  different 
temperatures  for  the  DMFC  operating  with  air  at  the  cathode  and 
at  a  methanol  feed  concentration  of  1.0  M,  while  power  density  is 
plotted  as  a  function  of  current  density  in  Fig.  3(b).  Although  the 
plots  are  qualitatively  similar  to  those  in  Fig.  2,  there  are  important 
differences.  First  of  all,  there  is  a  steady  increase  with  tempera¬ 
ture  in  cell  voltage  as  well  as  power  density  for  a  given  current 
density  throughout  the  range  of  current  densities.  This  is  because 
the  higher  temperatures  significantly  enhance  electrode  kinetics, 
as  well  as  transport  processes,  both  proton  transport  in  the  PEM 
[5],  as  well  as  diffusion  across  the  GDLs.  The  increase  in  cell  voltage 
and  power  density  with  temperature  at  lower  current  densities  is 
a  clear  indication  of  improving  electrode  kinetics  at  both  the  anode 
and  the  cathode.  The  reduction  of  slope  in  the  polarization  plot  in 
Fig.  3(a)  at  the  intermediate  current  densities  indicates  a  reduc¬ 
tion  in  the  Ohmic  resistance  of  the  PEM.  Finally,  the  increase  in 
limiting  current  densities  with  temperature  is  a  clear  indication  of 
improvement  in  the  transport  across  the  GDL.  In  short,  the  higher 
the  operating  temperature  in  the  range  investigated,  the  better, 


Table  2 

Experimental  conditions  of  Chiu  et  al.  [4,20]. 


Item 

Description 

Commercial  MEAs 

Anode  catalyst 

Cathode  catalyst 

Pt  loading 

Graphite  channel 

Anode  flow  rate  (liquid  feed),  VA 
Cathode  flow  rate  (Air),  Vc 
Operating  temperatures,  T  (I<) 
Methanol  feed  concentrations, 

Qvie.in  (M) 

35  mm  x  35  mm  duPont,  with  Nafion®  117 

40  wt%  Pt-Ru/C,  Pt/Ru  wt  ratio  1:1 

40wt%  Pt  1C 

2  mg  cm-2  at  both  anode  and  cathode 

Single  serpentine,  1  mm  wide  x  1  mm  deep 

5  ml  min-1 

150  ml  min-1 

303,  308,  313,  318,  323,  328,  333 

0.75,  0.875, 1.0, 1.125, 1.25, 1.375, 1.5,  2.0 

provided  it  is  not  so  high  that  the  cathode  is  asphyxiated  by  the 
rapidly  rising  vapor  pressure  of  water. 

Further,  methanol  crossover  escalates  with  temperature  as  well, 
so  that  we  might  expect  that  there  is  an  optimum  temperature  for 
a  given  methanol  feed  concentration. 

The  maximum  cell  power  density,  Pm ax,  as  a  function  of  tem¬ 
perature  and  feed  composition,  as  shown  in  Figs.  2(b)  and  3(b)  can, 
in  fact,  be  determined  analytically  by  taking  the  derivative  of  the 
power  density  (Eq.  (37))  with  respect  to  current  density,  and  setting 
it  equal  to  zero  to  determine  the  corresponding  optimum  current 
density,  which,  when  used  back  in  Eq.  (37),  provides  Pmax .  Such 
information  is  crucial  for  developing  control  strategies  for  a  DMFC 
power  plant,  e.g.,  for  real-time  control  of  DMFC. 

3.2.  Crossover  and  its  impact  on  electrode  overpotentials  and  OCV 

Methanol  crossover  is  the  central  issue  in  DMFCs  that  is  known 
to  not  only  represent  fuel  loss  but  also  increase  the  cathode  overpo- 
tential.  To  investigate  methanol  crossover  further,  Fig.  4  provides  a 
comparison  of  the  experimental  results  of  Chiu  et  al.  with  the  pre¬ 
dicted  methanol  crossover  current  density  (Eq.  (27))  as  a  function  of 
the  fuel  cell  current  density  at  different  operating  temperatures  for 
a  1 .0  M  methanol  feed  and  using  the  set  of  parameters  in  Table  1 .  It 


i  (A/cm2) 


Fig.  4.  Comparison  of  theoretical  predictions  with  experiments  of  Chiu  et  al.  for 
methanol  crossover  current  density  versus  current  density  at  different  tempera¬ 
tures  for  a  DMFC  operating  withl.O  M  methanol  feed  and  air  at  the  cathode.  Table  2 
provides  the  other  experimental  conditions. 
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Fig.  5.  Comparison  of  theoretical  predictions  with  experiments  of  Chiu  [4]  for 
methanol  crossover  current  density  versus  current  density  at  different  methanol 
feed  concentrations  for  a  DMFC  operating  at  323  K  with  air  at  the  cathode.  Table  2 
provides  the  other  experimental  conditions. 


V(V) 


Fig.  6.  Predicted  anode  and  cathode  overpotentials,  including  contributions  of 
methanol  crossover,  of  a  DMFC  operating  at  50 °C  with  a  2M  methanol  feed  and 
air  at  the  cathode. 


is  noteworthy  that  crossover  current  density  is  of  the  same  order  as 
fuel  cell  current  density  even  at  this  low  methanol  feed  concentra¬ 
tion.  Methanol  crossover  is,  of  course,  contributed  by  both  diffusion, 
dependent  on  methanol  concentration,  and  electro-osmotic  drag, 
dependent  on  current  density. 

At  zero  current  density,  or  under  OCV  conditions,  the  methanol 
crossover  current  density  is  the  highest,  corresponding  to  Eq.  (28), 
and  accounts  for  the  significant  drop  of  OCV  below  the  thermody¬ 
namic  voltage  of  DMFC.  It  is  seen  in  Fig.  4  that  crossover  current 
density  at  OCV  increases  sharply  with  temperature.  As  the  DMFC 
current  density  increases,  however,  more  of  the  methanol  is  con¬ 
sumed  within  the  ACL  and,  thus,  less  is  available  for  diffusion 
across  the  PEM,  so  that  the  crossover  current  density  declines 
(Fig.  4),  eventually  approaching  zero  as  the  fuel  cell  current  den¬ 
sity  approaches  the  limiting  current  density.  Nonetheless,  there  is 
a  sharp  increase  in  the  methanol  crossover  at  all  current  densities 
with  temperature. 

There  is  also  a  qualitatively  similar  relationship  between 
methanol  crossover  current  density  versus  fuel  cell  current  den¬ 
sity  for  varying  methanol  feed  concentrations  at  a  given  operating 
temperature,  as  shown  in  Fig.  5.  It  is  seen,  once  again,  that  the 
model  does  a  nice  job  of  simulating  the  experimental  data  of  Chiu 
[4]  without  any  additional  fitting. 

Since  the  model,  thus,  appears  robust,  let  us  use  it  to  further 
explore  the  effect  of  crossover  on  electrode  overpotential  as  well 
as  on  OCV.  Thus,  the  model  is  used,  along  with  the  parameters  in 
Table  1 ,  to  theoretically  predict  the  cathode  as  well  as  anode  polar¬ 
ization  as  shown  in  Fig.  6,  for  the  experimental  conditions  of  Chiu 
et  al.  [20]  (Table  2),  and  for  an  assumed  methanol  feed  concentra¬ 
tion  of  2  M  and  an  operating  temperature  of  323  K. 

The  upper  curve  in  Fig.  6  represents  the  cathode  polarization 
because  of  kinetic  and  diffusion  limitations  of  ORR,  r]C  =  i)c,i<  +  ilc,D> 
while  the  yellow  sliver  below  it  is  the  contribution  of  methanol 
crossover  to  the  cathode  overpotential,  r]CtMe,x-  Thus,  the  upper 
shaded  area  represents  cathode  power  consumption  for  a  H2-02 
fuel  cell,  to  which  when  is  added  the  yellow  area  representing  the 
polarization  and  corresponding  power  consumption  because  of 
methanol  crossover,  the  total  cathode  polarization  loss  due  to  a 
DMFC  cathode  results.  It  is  seen  that  while  the  ORR  sluggishness 
is  the  dominant  cause  of  cathode  polarization,  the  contribution 
of  the  methanol  crossover  is  significant  especially  at  the  lower 
current  densities,  in  accord  with  experimental  observations 
[43].  Further,  when  the  ORR  is  diffusion  limited  at  high  current 


densities,  although  not  shown  in  Fig.  6,  again  the  contribution  of 
methanol  crossover  to  cathode  overpotential  can  become  signif¬ 
icant,  as  it  reduces  the  effective  cathode  limiting  current  density 
because  of  oxygen  consumption  by  the  methanol  that  has  crossed 
over. 

The  lower  region  in  Fig.  6  represents  the  anode  polarization 
losses.  Again  the  yellow  sliver  in  the  figure  represents  the  effect 
of  methanol  crossover.  Thus,  there  is  not  an  insignificant  anode 
polarization  loss  due  to  methanol  crossover.  The  orange  region 
below  it  represents  kinetic  and  diffusion  limitations  of  the  MOR, 
rlA  =  rlA,K  +  rlA,Di  determined  by  setting  the  term  representing  the 
crossover  current  density  of  methanol  equal  to  zero  in  the  anode 
overpotential,  Eq.  (30). 

It  is,  thus,  noteworthy  that  the  methanol  crossover  is  predicted 
to  increases  the  polarization  of  the  anode  as  well,  in  addition,  of 
course,  to  that  of  the  cathode,  something  that  has  not  apparently 
been  discussed  in  the  literature  so  far.  The  limiting  current  den¬ 
sity  of  the  cell  in  the  example  in  Fig.  6  is  determined  by  iAL,  which 
in  turn  is  determined  mainly  by  the  methanol  feed  concentration 
[44].  Finally,  the  colorless  region  between  the  cathode  and  the 
anode  polarization  regions  in  Fig.  6  represents  the  useful  electric 
power  produced  in  the  DMFC,  the  vertical  distance  between  the 
two  overpotential  curves  being  the  observed  voltage,  V. 

Another  potential  issue  affecting  DMFC  performance  that  has 
not  so  far  been  discussed  in  the  literature  is  the  anode  polariza¬ 
tion  because  of  oxygen  crossover,  illustrated  in  Fig.  7.  Although  in 
a  well-assembled,  air-fed  DMFC,  the  02  crossover  current  is  small 
(<1  mAcnrr2),  it  could  be  non-negligible  when  the  anode  chamber 
is  open  to  air,  as  in  a  passive  DMFC  to  vent  C02,  or  when  account  is 
taken  of  the  02  dissolved  in  an  air-equilibrated  methanol  feed.  As 
Fig.  7  shows,  ingress  of  02  at  the  anode  can  not  only  substantially 
impact  the  OCV,  but  also  the  entire  polarization  plot  as  well  as  the 
diffusion  restricted  anode  limiting  current  density. 

According  to  Eq.  (38),  the  thermodynamic  voltage  in  a  DMFC 
under  typical  operating  conditions  is  about  1.18  V.  However,  the 
actual  OCV  is  only  about  half  as  much  [20],  depending  upon 
the  operating  temperature  and  the  methanol  feed  concentration. 
This  is  primarily  because  of  the  combined  effect  of  the  methanol 
crossover,  and  the  sluggish  ORR  in  oxidizing  the  methanol  at  the 
cathode  (Fig.  6).  As  the  methanol  concentration  increases,  methanol 
crossover  becomes  more  significant  (Fig.  5),  reducing  the  OCV  even 
further. 
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Fig.  7.  The  predicted  effect  of  oxygen  crossover  current  density  on  anode  overpo¬ 
tential. 
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Fig.  9.  A  comparison  of  predicted  heat  flux  dissipation  and  electric  power  density 
produced  versus  current  density  in  a  DMFC  operating  at  323  K  with  a  2  M  methanol 
feed  and  air  at  the  cathode. 


The  OCV  can,  in  fact,  be  predicted  simply  by  setting  i  =  0  in  the 
DMFC  polarization  expression,  Eq.  (38),  resulting  in 


OCV  =  V0  - 


RT 

UAVA,e-F 


sinh  1 


1  (  *X,O2/*A,0  A  f,  *X,M e,l\  1 

2  -ix,02/^A,L  )  \  (a, I  )  J 


+  - 


RT 


acVc,e-F 


sinh 


(x,Me,o/  ic,0 
1  -  iX,Me,0 /k,L 


(59) 


i.e.,  both  oxygen  and  methanol  crossover  can  affect  the  OCV.  Thus, 
the  OCV  is  plotted  in  Fig.  8  as  a  function  of  temperature  as  well 
as  feed  methanol  concentration  for  a  DMFC  operating  on  air  at  the 
cathode.  Clearly,  the  OCV  declines  monotonically  with  increasing 
methanol  concentrations  at  all  temperatures. 

The  dominant  effect  on  OCV  is,  of  course,  that  of  the  methanol 
crossover  under  zero  current  conditions,  i.e.,  i’x.Me.o*  as  given  by 
Eq.  (28).  From  Eqs.  (28)  and  (59),  thus,  for  the  case  when  the  oxygen 


Fig.  8.  Predicted  OCV  as  a  function  of  operating  temperature  (300-350  K)  and  feed 
methanol  concentration  (0.1 -5.0  M,  or  0.1  x  10-3  to  5.0  x  10-3  mol  cm-3)  for  an  air 
fed  DMFC. 


crossover  may  be  neglected,  the  OCV  may  be  written  as 


OCV  =  Vo  +  ^ - =sinh-1 

acvc,e-F 


_ ix,Me,l/2ic,0 _ 

1  +  ix,Me,l(( l/*A,r)  -  (l/ic,l)) 


(60) 


While  the  OCV  declines  monotonically  with  concentration  at  a 
given  temperature  (Fig.  8),  it  is  seen  that  at  higher  methanol 
feed  concentrations,  the  OCV  is  not  a  monotonic  function  of 
temperature,  although  that  is  the  case  at  the  lower  methanol  con¬ 
centrations.  The  electrode  kinetics  at  the  anode  as  well  as  the 
cathode  improve  with  temperature,  reducing  the  electrode  over¬ 
potentials.  On  the  other  hand,  methanol  diffusion  through  the 
membrane  increases  with  temperature,  as  does  the  crossover  with 
temperature.  Even  more  significantly,  the  vapor  pressure  of  water 
increases  quickly  with  temperature  (Table  1 ),  smothering  the  cath¬ 
ode  and  consequently  reducing  the  cathodic  diffusion  restricted 
limiting  current  density,  iCL. 


3.3.  DMFC  Efficiency 


The  enthalpy  of  combustion  of  the  fuel,  AH,  that  is  not  con¬ 
verted  into  useful  electric  power  (Eq.  (37))  because  of  the  various 
irreversibilities  in  the  anode,  cathode,  and  the  PEM,  is  dissipated 
by  the  DMFC  as  low-grade  heat  [45],  which  must,  consequently,  be 
managed  for  maintaining  the  desired  operating  temperature.  The 
heat  dissipation  flux  (W  cm-2  MEA)  is  determined  from 

Q  =  Vmax(i  +  (x,Me  +  *X,02  )  —  P  (61 ) 

where  P  is  given  by  Eq.  (37),  and  Vmax  =  -A H/(ne-F).  For  DMFC, 
V£ax  =  1-255V  under  standard  conditions,  while,  of  course,  the 
thermodynamic  potential  V°  =  -  A G/ne-F  =  1.214  V  under  stan¬ 
dard  conditions,  i.e.,  the  thermodynamic  efficiency  under  standard 
conditions,  £°  =  V°/V°ax  =  1.214/1.255  =  96.7%. 

The  heat  flux  generated  versus  current  density  is  plotted  in 
Fig.  9  for  illustrative  purposes  for  the  operating  conditions  of 
Fig.  6.  It  is  obvious  that  a  substantial  amount  of  heat  is  generated 
in  a  PDMFC,  while  only  a  relatively  small  fraction  of  the  energy  of 
the  fuel  is  converted  into  electric  power,  much  as  Eccarius  et  al. 
[46]  have  reported.  Thus,  the  DMFC  suffers  from  a  low  efficiency, 
clearly  because  of  the  large  electrode  overpotentials  coupled  with 
methanol  crossover,  which  reduces  the  current  efficiency.  The 
generated  heat  can,  in  fact,  be  used  to  raise  the  operating  temper¬ 
ature  of  a  passive  DMFC  [47]  and,  thus,  improve  performance,  as 
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Fig.  10.  Comparison  of  theoretical  predictions  of  efficiency  versus  current  density, 
with  experimental  efficiency  calculated  from  the  data  of  Chiu  et  al.  [20]  at  differ¬ 
ent  temperatures  for  a  DMFC  operating  with  a  1.5  M  methanol  feed  and  air  at  the 
cathode.  Table  2  provides  the  other  experimental  conditions. 


Fig.  11.  Comparison  of  theoretically  predicted  petal-shaped  plots  of  efficiency 
versus  power  density,  with  experimental  efficiency  calculated  from  the  data  of  Chiu 
et  al.  [20]  at  different  temperatures  for  a  1 .5  M  methanol  feed  and  air  at  the  cathode. 
Table  2  provides  the  other  experimental  conditions. 


shown  by  Casalegno  et  al.  [48].  However,  as  discussed  above,  too 
high  a  temperature  also  enhances  methanol  crossover,  and  further 
stifles  the  cathode  because  of  the  higher  vapor  pressure  of  water 
within  the  cathode  chamber. 

The  overall  DMFC  efficiency  [49,50]  may  be  determined  from 
a  product  of  the  thermodynamic,  current,  and  voltage  efficiencies, 
i.e., 

e,Wy,(  .  ■  .  )(£) 

\Vmax/  y  *  +  lx,Me  +  lX,02  )  V  ^0  / 


Vmax(i  +  *X,Me  +  *X,02  ) 

Fig.  10  shows  a  comparison  of  theoretical  predictions  with  experi¬ 
mental  efficiencies  at  different  temperatures  for  a  1.5  M  methanol 
feed  and  air  on  the  cathode,  calculated  based  on  Eq.  (62)  and  uti¬ 
lizing  the  experimental  data  of  Chiu  et  al.  [20]  for  power  density, 
P,  current  density,  i,  and  the  crossover  current  density,  ix,Me»  while 
assuming  that  the  oxygen  crossover  current  density,  ix, o2  0.  A 
similar  comparison  is  provided  in  Fig.  1 1  for  DMFC  efficiency  versus 
power  density  at  different  temperatures  for  a  1.5  M  methanol  feed 
and  air  on  the  cathode. 

The  maximum  efficiency  observed  for  the  DMFC  conditions 
investigated  is  around  24%,  which  does  not  appear  to  vary  much 
with  temperature,  although  the  peaks  becomes  broader  with  tem¬ 
perature,  i.e.,  higher  efficiencies  can  be  obtained  over  a  broader 
range  of  operating  current  densities  with  temperature.  This  is  use¬ 
ful  in  applications  where  the  load  can  vary  significantly. 

The  peaks  in  efficiencies  observed  in  Fig.  11  as  a  function  of 
power  density  appear  even  sharper  in  the  petal-shaped  plots.  The 
reason  these  petal  plots  go  through  the  origin  is  that  the  efficiency 
6  =  0  when  P=  0,  both  under  OCV  conditions,  when  i  =  0,  as  well 
as  under  limiting  current  density  conditions,  when  V=0.  So,  as 
one  traverses  from  OCV  to  limiting  current  density  (clockwise  on 
the  petal),  the  efficiency  first  increases,  peaks,  and  then  declines, 
returning  back  to  zero.  The  agreement  with  theory  is  once  again 
quite  good.  It  is  seen  that  there  is  only  a  narrow  range  of  opera¬ 
tion  indicated  if  peak  power  density  as  well  as  peak  efficiency  are 
desired.  Such  plots  and  relationships  can,  thus,  be  very  useful  for 
algorithmic  control  of  DMFCs  to  operate  in  the  region  of  peak  effi¬ 
ciencies  and  power  densities,  which  is  clearly  important  in  portable 
applications. 


Finally,  although  not  discussed  in  this  paper,  the  model  pre¬ 
sented  here  is  useful  in  exploring  design  aspects  of  DMFCs,  e.g., 
effect  of  the  catalyst  loading  and  particle  size,  and  membrane 
thickness  [51],  etc.  In  general,  thicker  membranes  perform  bet¬ 
ter  at  lower  current  densities  and  generate  higher  OCV  due  to 
lower  crossover,  although  thinner  membranes  exhibit  better  per¬ 
formance  at  higher  current  densities  because  of  smaller  electrolyte 
resistance.  Of  course,  the  model  may  be  used  for  diagnostic  pur¬ 
poses  as  well,  e.g.,  a  decline  in  OCV  may  be  an  indication  of 
membrane  degradation/thinning.  The  interfacial  resistance  in  the 
above  analysis  has  been  assumed  to  be  zero.  While  this  may  be  valid 
for  a  well-fabricated  MEA,  an  increase  in  the  resistance  with  time 
might,  for  instance,  indicate  delamination  of  the  layers,  or  simply 
a  poorly  fabricated  MEA. 

4.  Conclusions 

The  direct  methanol  fuel  cell  is  very  attractive  for  portable  and 
backup  power,  and  is  starting  to  make  inroads  into  some  of  these 
markets,  e.g.,  as  an  on-the-go  forklift  battery  charger.  However,  its 
high  cost  and  relatively  unimpressive  power  density  remain  sig¬ 
nificant  barriers.  Operating  DMFC  under  optimal  power  density 
and  efficiency  requires  a  clear  understanding  of  how  the  var¬ 
ious  design  and  operational  parameters  affect  its  performance. 
We  have,  thus,  developed  a  transparent  albeit  rigorous  analytical 
model  based  on  some  key  judicious  but  defensible  assumptions. 
The  resulting  analytical  expressions  for  crossover  current  density, 
DMFC  polarization,  and  power  density  are  in  terms  of  quantities 
with  discernable  physical  meaning,  so  that  their  effect  on  DMFC 
performance  can  be  better  appreciated. 

The  model,  with  mainly  two  fitted  parameters,  has  been  rigor¬ 
ously  tested  against  an  extensive  and  consistent  data  set  on  DMFC 
performance,  and  found  to  be  in  broad  agreement.  Thus,  the  model 
accurately  predicts  the  influence  of  methanol  crossover  on  OCV 
as  well  as  on  the  overall  cell  performance  under  a  variety  of  feed 
concentrations  and  temperatures.  The  fuel  cell  power  density  and 
efficiency  under  different  operating  can  be  accurately  predicted  as 
well.  Consequently,  the  model  can  be  used  for  developing  control 
algorithms  for  operating  under  optimal  conditions,  e.g.,  by  carefully 
controlling  online  the  methanol  feed  concentration.  It  can  also  be 
used  for  diagnostics.  Finally,  the  model  should  help  in  improved 
understanding  of  the  DMFC. 
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